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SUMMARY 


Numerical models for two-dimensional turbulent boundary 
layers over wavy surfaces have been investigated. Computations 
for wavy v/all boundary layers indicate possibilities of overall 
drag reduction in a parameter range of the geometry of the wall. 
The correction technique using integral methods for analyzing 
arbitrary surfaces has been found to be unsuitable for some cases 
of interest in drag reduction: a Navier-Stokes solver for v;avy 

walls has, hence, been built to test these problems. Test 
results of the Navier-Stokes solver indicate that the solution 
techniques are accurate enough to handle complex geometries and 
steep variations in fluid properties. 



LIST OF SYMBOLS 


a amplitude of the wavy surface 

K curvature of the wavy surface 

k turbulent kinetic energy 

p pressure 

R,r Radius of curvature of the wavy surface 
w 

t time 

u velocity in the longitudinal direction 

V velocity in the normal direction 

X distance along the longitudinal direction 
y normal distance 

6 boundary layer thickness 

■k 

6 momentum deficit in the boundary layer 

6 . . Kronecker delta 
ID 

X wave length 

dynamic viscosity in laminar and turbulent cases 

' ,a ' aoparent normal stresses due to turbulent flow 
X y 

T apparent shear stress due to turbulent flow 

2 

V Laplacian 
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INTRODUCTION 


The energy crisis of recent years has brought to focus the 
need for testing both conventional and unconventional ideas of 
turbulence control. Some candidate concepts of this type are 
being studied at NASA Langley Research Center and include compli- 
ant surfaces, wavy geometries and turbulent eddy breakup devices. 
The present contract work was initiated to provide logistical 
support for the research in turbulence control and drag reduction 
at the Fluid Mechanics Branch, LaRC, NASA. 

The current state of experimental and theoretical research 
on wavy geometries has been reviewed in ref. 1. Briefly stated 
there are two condidate experiments (i.e. references 2 and 3) that 
shov; substantial skin friction drag reduction in flows over wavy 
surfaces. However the total drag for these surfaces is above the 
flat plate value due to the enormous pressure drag contribution. 

In reference 1 details of a procedure for predicting skin 
friction drag and the pressure drag over v/avy surfaces, have been 
presented. The present research effort, based on this work was 
directed to develop algorithms for the design of v;avy wall experi- 
ments and to predict the performance of these geometries. Details 
of this effort are docum.ented in Section 2. 

In section 3 details of the effort to develop a two-dimension- 
al Navier-Stokes solver for wavy geometries are presented. In 
section 4 a review is presented of the feasibility study of candi- 
date concepts and directions for future research. 
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2 . 


FLOW OVER WAVY SURFACES 


turbule: 


Turbulent flows over wavy geometries undergo two major 
effects which are absent in flat geometries, (a) oscillatory 
curvature effects, and (b) oscillatory pressure gradient effects. 
For waves of small curvature in thick turbulent boundary layer 
the v/avy surface may be viewed as a perturbation from flat 
geometry and estimates on flov; modifications are sufficiently 
easy to predict. However v/here certain parameters (^) (i.e. 

wave number) and (6/R) (i.e. nondimensional curvature) are moder- 

ate the problem can no longer be viev/ed as a small perturbation 
of the flat boundary layer case. The mean boundary layer equations 
for two-dimensional incompressible turbulent flows of non- 
negligible surface curvature in body-fitted coordinates are given 
below : 

Continuity : 


k P “ / ^ (p V + p- V )| = 0 

(la) 


Mean Momentum Equations : 
1 


(1 + Ky) 


p U + (p V + P ' V ») + (P V H P » V »)* 

o X ^ Y 


uK 


9P + ^ + 2K 


1 + Ky 3x 3y (1 + Ky) 


(Ky + 1) 

T 


K 


1 + Ky 3y 
K 

1 + Ky 


9 , , 

^ (y u) 


- u u 


K 


U + Kyj 


P u 


= 

3y 


(Ib) 


( Ic) 



The presence of equation Ic is due to non-negligible 
gradients of pressure introduced by the curvature. The conven- 
tional boundary layer assumption (-^ =0) is thus inappropriate. 

The equation system I given above is difficult to solve due to 
the presence of equation Ic. The closure model for turbulence 
must include the non-equilibrium nature of the flow over the 
waves and their relaxation to equilibrium when the curvature 
and pressure gradient effects are applied to the flovz. 

In reference 1 a procedure for solving the equations of 
motion for the turbulent flov; over wavy surfaces was presented. 

The procedure adopted there is to use an existing two-dimensional 
boundary layer code (i.e. the Beckwith-Eushnell compressible 
boundary layer program of ref. 4) with a specified input pressure 
field and suitable modifications to length scales due to curvature 
and pressure gradient. The prescription of P(y) means that to 
first order the vertical momentum equation can be ignored. Thus 
one is left with solving a boundary layer program to obtain the 
required results such as skin friction, velocity profiles, etc. 

The approach to wavy wall turbulent flows therefore involves: 

(i) evaluation of pressure distribution for input to the 
program, and 

(ii) use of the modified Beckwith-Bushnell program for obtain- 
ing the required solution. 

The pressure distribution over the wavy surface is evaluated 
using an inviscid perturbation method of the Miles type discussed 
in ref. 5. An algorithm which is third order accurate in space 
has been developed to obtain the pressure distribution for this 
case. Unfortunately, this algorithm is of the inviscid type and 
hence no information of the pressure phase shift can be obtained 
by this method. The phase shift of the pressure is currently 
evaluated using empirical data as shov7n in fig. 1. 
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A procedure for designing wavy wall surfaces with total drag 
reduction is given in ref, 1. However, extensive testing of thi 
approach on moderately curved wavy surfaces indicates that the 6 
correction approach does not worlc for these problems. The failure 
of this first order approach suggests that the correct way to 
approach the problem should be by solving the full Navier-Stokes 
equations on the computer. In the next section details of such 
a procedure are described. 
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3. 


NAVIER-STOKES ALGORITHM FOR WAVY SURFACES 


The governing equations for a two-dimensional flow over a 
wavy surface in a Cartesian x-y coordinate system are: 


3u 

3x 


+ 


3v 

3y 


= 0 


3u , 3u 


+ V 


3u 

3y 


-1 

P 


lE 

3x 


+ V 


L^7 


+ 


3^u"^ 

3x2 u 




3a 

X 

3x‘ 


+ — ) 
3y ^ 


3v dy_ _ ^ 

3t ^ 3x ^ 9y P 3y 


+ V 


1-0x2 3y2 


] 


, ^ * do 

-H i ( ^ 

p ^ 3x 3y 


(II) 


a a 

The modelling of viscous terms due to turbulence, ( x, t, y ) 
is done by using eddy vicosity (i.e. through ad hoc length scales 
v/hich seem to work very v/ell for moderate [ja/X) , (6/R) case^ ) , 

obtained from the 2D modified Beckwith-Bushnell code (ref. 4). 
Thus , 


X 


= -p u 


= -p u ' V 


a E -p V 

y 


.2 


and - 


(III) 


p Uj = ( 


3u . 


3u 


t ' X . 3 X . 

11 


) - f Pk 6,. 
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The relation for k is 


where is turbulent eddy viscosity, 
provided by empirical data as: 


u’ 


v' 

k 


0.15 


(IV) 


Since is a function of both x and y, equation system II has 
to be solved very carefully to obtain the required accuracy. 


Proposed Method 

The method for solving the Navier-Stokes equations is by 
converting the semi-infinite domain of analysis into a rectangular 
box by using a conformal mapping program which uses FFTs and is 
hence very inexpensive. In fig. 2 the mapping is shown schemat- 
ically. The contra and covariant derivatives connecting the 
derivatives in the mapped plane and in the untransformed plane 
are, evaluated very cheaply and stored in a preprocessing routine. 
The dependent variables are expanded in Fourier-Chebychef f modes. 
The solution technique used guarantees spectral accuracies in 
space thereby limiting the core storage to some manageable level. 
Velocity splitting is used. A schematic of the computation is 
shown in fig. 3. A Richardson extrapolation technique is used to 
obtain second order accurate solution in tirr.e. The spatial accura- 
cy is N fold where N is the number of grids used for solution. At 
present the analysis is confined to over o.nly one wavelength to 
keep the expenses of running the code minimal. 

Initial test runs using the above spectral code indicate 
spectral accuracy and that the method is second order accurate in 
timie . 
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4, DISCUSSION OF RESULTS 


Wavy Wall Drag Reduction 

Figures 4 and 5 are the cumulative findings of a large number 
of test runs using the modified Beckv/ith-Bushnell program. The 
pressure phase shift is input by relying on fig. 1; the pressure 
decay by calculations using the algorithms discussed in Section 2. 
The results of fig, 4 indicate that (a/A) of 0.004 gives the 
largest reduction in drag. Figure 5 indicates that again for this 
waveshape the equilibrium value of total drag is attained after 
about 40 wavelengths and the net drag reduction is not improved 
any further farther downs treacA, 

Navier-Stokes Test Results 

The Navier-Stokes code that has been developed is in nodular 
form. The various modules and their interconnection are shov/n 
in figure 6. The Navier-Stokes algorithm was initially tested to 
check for any bugs in the various models. Two of the important 
nodules that needed to be tested thoroughly were the Poisson 
solver and viscous stabilization solvers. Below the results of 
these tests are presented. 

Test No. 1 : 

Stokes Equation 

The equation system given below v/as solved: 


u^ = 

-p 

+ u 

+ u 


u 

X 

XX 

yy 


v. = 

-p 

+ V 

+ V 

(V) 

t 

u 

XX 

yy 

u + 

V 

= 0 




y 





vjith u = sin 2 x e " 

V = + 2 Cos 2 X e as the initial velocity field. (VI) 
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The time evolution indicated decay of the velocity profile as 
-3t 

e and the scheme was found to be second order accurate in 
time. The code was run with meshes of 6x8, 6x16, 32x32 and 
64x84 and also Sxl6, 8x32 and 3x64. It v;as found that if the 
Fourier content involves only the first four wavs numbers only 
8 points are required in the x directions and v/ith 32 and 64 
modes in the Tchebycheff direction snectral accuracies can be 
obtained. This test qives confidence that the heart of the 
Mavier-Stokes solver, i.e. the Poisson solver, is v;orking very 
we 1 1 . 

Test No. 2: 

u^ = g V u 

vjith u = sin X e (VII) 

3 1 * 

and u (x,o) = e sin x 

Again the test results indicate that the viscous updates 
are solved exactly and that the viscous stabilization module is 
working well. 
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CONCLUSIONS 


A Wavier-Stokes solver for flov/ over two-dimensional wavy 
surfaces has been constructed. The Navier-Stokes solver is 
extensively debugged and can be used for solving flow over 
arbitrarily wavy surfaces. The Navier-Stokes code can use time 
steps of about four orders of magnitudes higher than the explicit 
CFL lirait/ it is second order accurate in time and infinite accu- 
rate in space. 

The computation of wavy flov/ over moderate amplitude waves 
indicates some total drag reduction in a parameter range pre- 
viously left untouched by experimenters which opens up exciting 
possibilities for future research. 

The author would like to thank Mr. A. Cary for his con- 
tribution to the work on wavy wall computations. The author 
v/ould also like to thank Prof. S. A. Orszag for providing 
guidance in setting up the Navier-Stokes code. 
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WAVY WALL 


Schematic of Mappings 


Figure 2 











Schematic of Calculation 
Procedure for Wavier Stokes Problem 


Figure. 3 
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AVERAGED OVER 50 WAVES 



Figure 4 
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WAVE NO. 

Analytical predictions of the average drag per wave 
for several shallow sine-wave walls. 

Figure 5 
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